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ABSTRACT 

We present a radio continuum study of the pulsar wind nebula (PWN) DA 495 (G65.7+1.2), including 
images of total intensity and linear polarization from 408 to 10550 MHz based on the Canadian Galactic 
Plane Survey and observations with the Effelsberg 100-m Radio Telescope. Removal of flux density 
contributions from a superimposed H II region and from compact extragalactic sources reveals a break 
in the spectrum of DA 495 at 1.3 GHz, with a spectral index a = —0.45 ± 0.20 below the break and 
a = —0.87 ± 0.10 above it {S^ ex v"). The spectral break is more than three times lower in frequency 
than the lowest break detected in any other PWN. The break in the spectrum is likely the result of 
synchrotron cooling, and DA 495, at an age of ^20,000 yr, may have evolved from an object similar to 
the Vela X nebula, with a similarly energetic pulsar. We find a magnetic field of ^1.3 mC inside the 
nebula. After correcting for the resulting high internal rotation measure, the magnetic field structure is 
quite simple, resembling the inner part of a dipolc field projected onto the plane of the sky, although 
a toroidal component is likely also present. The dipole field axis, which should be parallel to the spin 
axis of the putative pulsar, lies at an angle of ~50° east of the North Celestial Pole and is pointing 
away from us towards the south-west. The upper limit for the radio surface brightness of any shell- type 
supernova remnant emission around DA 495 is Yji ghz ^ 5.4 x 10~^^ Watt m~^ Hz"-'^ sr"-'^ (assuming a 
radio spectral index of a = —0.5), lower than the faintest shell-type remnant known to date. 

Subject headings: ISM: individual (G65.7-I-1.2), magnetic fields, polarization, supernova remnants 



1. INTRODUCTION 

The Crab Nebula, the remnant of the supernova explo- 
sion of 1054, is an astrophysical laboratory of immense 
value where we can study the development of a young pul- 
sar and the synchrotron nebula generated by the particles 
that it injects. The Crab Nebula has become the archetype 
of a class of objects earlier known as filled-center, plerionic, 
or Crab- like supernova remnants (SNRs), but now called 
pulsar wind nebulae (PWNe) in recognition of their energy 
source. PWNe are a minority in SNR catalogs, less than 
10% in number, but if we include the composite SNRs, syn- 
chrotron nebulae centrally placed within the more usual 
SNR shell, their ranks sweU to 15% of the Galactic SNR 
catalog (Green 2004). 

PWNe are characterized by a flat radio spectral index, 
a, in the range > a > —0.3 (where the flux density S at 
frequency v varies as S cc v"') and a central concentration 
of emission gradually declining to the outer edge. Their 
appearance reflects the central energy source and their ra- 
dio spectrum reflects the energy spectrum of the injected 
particles. The more common shell SNRs exhibit steep 
outer edges and steep radio spectra (generally a < —0.3). 

We have historical records of fewer than ten supernova 
events (Stephenson & Green 2002). By chance two of those 
have left PWNe, the Crab Nebula and 3C 58, with no de- 
tectable trace of accompanying shells, and one a compos- 
ite remnant, Gil. 2— 0.3. Thus we have three examples in 
which we can closely follow the evolution of young PWNe, 
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but the relatively small number of known PWNe has pro- 
vided little observational basis for our understanding of 
the evolution of these objects through later stages. This 
paper presents and analyzes new observations of DA 495 
(G65.7+1.2), which seems likely to be a PWN of advanced 
age. 

DA 495, while an outlier among PWNe, is without doubt 
a member of that class. Its r adio appearance, first seen i n 
the high-resolution images of lLandecker fc Caswelll (|1983[ ). 
has the signature of a PWN, a central concentration of 
emission smoothly declining with increasing radius, with 
no evidence of a surrounding shell (although there is a de- 
pression in intensity near its center). The only indication 
that it might not be a PWN is its mea sured radio spectral 
index, a ~ —0.5 (jKothes et al.ll2006bl ). much steeper than 
is usual for this class of objects, and more closely resem- 
bling that of a shell-type SNR. Nevertheless, shell rem- 
nants have steep outer edges where kinetic energy is being 
converted to radio emission near the shock front, while 
in DA 495 the emissivity distribut ion is centrally concen- 
trated (jLandecker fc Caswellll 19831 Fig. 3) indicating that 
the energy source is central. Neve rtheless, on the basis o f 
the central depression in intensity, IVelusamv et al.l (|1989| ) 
argue that DA 495 might be a shell remnant with an excep- 
tionally thick shell formed by a reverse shock interacting 
with an unusual amount of slowly moving ejecta. 

Recent analysis of X-ray data for DA 495, obtained from 
the ROSAT and ASCA archives, dramatically confirms the 
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PWN interpretation of this object by revealing a compact 
central obj ect surrounded by an exten ded non-thermal X- 
ray region (jArzoumanian ct al."2004'). Analysis of Chan- 
dra data (jArzoumanian ct al. 2008 ) shows a non-thermal 
X-ray nebula of extent ~30" surrounding the central com- 
pact object. The compact source is presumably the pul- 
sar (required to power the synchrotron nebula), and we 
shall refer to it as such (even though pulsed emission has 
not been detected). In this paper we present new radio 
images of DA 495 which include linear polarization infor- 
mation. Analysis of high-resolution radio images allows 
us to remove superimposed sources. After this the inte- 
grated spectrum of the source exhibits a spectral break 
characteristic of a synchrotron nebula, and the observed 
parameters allow an estimate of the strength of the mag- 
netic field within the nebula. Polarization images reveal 
the structure of the field within this PWN. 

2. OBSERVATIONS AND DATA ANALYSIS 

The DRAG observations at 408 and 1420 MHz are 
part of the Cana dian Galactic Plane Survey (CGPS, 
iTavlor et al.l 120031 ). Data from single-antenna telescopes 
were incorporated with the Synthesis Telescope data to 
give complete coverage of all structures down to the res- 
olution limit (det ails of the process can be found in 
I Taylor et al.l I2003D . The one exception is the polariza- 
tion data at 1420 MHz, for which no single-antenna data 
were available. However, the telescope is sensitive to all 
structures from the resolution limit ('^-^l') up to 45'. The 
absence of single-antenna data is not a concern for obser- 
vations of DA 495, whose fuU extent is < 30'. The CGPS 
polarization data were obtained in four bands (1407.0, 
1413.0, 1427.6, and 1434.5 MHz) each of width 7.5 MHz, 
allowing derivation of rotation measures (RMs) between 
those frequencies. 

The Effelsberg observations were carried out well be- 
fore this investigation began, and were retrieved from 
archived data. Observations at 4850 and 10550 MHz 
were made at different times, but in each case in excel- 
lent weather. The receiver systems, installed at the sec- 
ondary focus of the telescope, received both hands of cir- 
cular polarization, with two feeds at 4850 MHz and four 
feeds at 10550 MHz. Stan dard data-reduction software 
based on the N0D2 package (jHaslam et al.|[l983 ) was em- 
ployed throughout. Data reduction relied on "software 
beam switchi ng" to reject atmosp heric emission above 
the telescope (jMorsi fc Reichlll986D . Here every feed is 
equipped with its own amplifier to allow parallel recording 
of data. For each pair of feeds a difference map is calcu- 
lated in which disturbing atmospheric effects are largely 
subtracted. These dual-beam maps are restored to the 
eq uivalent single-beam m ap with the algorithm described 
bv lEmerson et al.l(|l979D . To reduce scanning effects, mul- 
tiple maps - coverages - of the field were observed at 
different parallactic angles. These scanning effects arise 
from instrumental offsets and drifts that produce different 
baselevels for each individual scan, a typical result from 
raster scanning. Individual coverages were combined us- 
ing the "plait" algorithm (jEmerson fc Grave I I1988D . by 
destriping the maps in the Fourier domain. We also used 
2695 MHz data fro m the Effelsberg Galactic Plane Survey 
(jReich et al.lll990r) . whose angular resolution is 4f3. 

Characteristics of the new observations of DA 495, made 



with the DRAG Synthesis Telescope and the Effelsberg 
100-m radio telescope, are presented in Table [TJ 

3. RESULTS 
3.1. DA 495 in total intensity 

Figure [T] presents total intensity images of DA 495 at 
408, 1420, 4850, and 10550 MHz as well as data at 
2695 MHz from the Effelsberg 11cm survey (j Reich et al.l 
Il990f ) and a 60^tm image from the IRAS Galaxy At- 
las (IGA, ICaoet al.| [T997). which is part of the CGPS 
database. The 1420 MHz image is completely con- 
sistent with the earlier DRAG image published by 
iLandecker fc Caswelll (|l983l) , although the new image has 
lower noise by a factor of '^ 3. All the radio images show a 
diffuse source of full extent about 25'. The high-resolution 
images at 1420, 4850, an d 10550 MHz show t h e cen - 
tral depression reported bv ILandecker fc Caswelll (|l983f ). 
Its position is RA(2000) = 19''52'"10.4^ DEC(2000) 
= -|-29°26f2, and it is now apparent that the pulsar 
(IWGA J1952.2-h2925) is not at the center of the d e- 
pression (as postulated bv ILandecker fc Caswelll (|1983[ )). 
but lies almost 2' east of it a t RA(2000) = 19^52 ' "17.0^ 
DEC(2000) = -|-29°25'53" (jArzoumanian et all I2008D . 
Clearly the explanation of the hole in the nebula cannot 
be quite as simple as a decline in emission from an aging 
pulsar. At larger radii, the emission drops off smoothly 
towards the outer edge without any sign of an outer shell 
at any frequency. 

At 408 MHz DA 495 is quite circular in appearance. 
However, an extension of DA 495 to the south-west is ev- 
ident at 1420 MHz and becomes increasingly prominent 
at higher frequencies. The IGA data shown in Figure [Tj 
reveal an infrared source coincident with this part of the 
radio source. The simplest interpretation is that an H II 
region is superimposed on the PWN, and in the analysis 
we assume this to be the case. 

3.2. The integrated spectrum of DA 495 

We derived the integrated emission of DA 495 from each 
of the radio images in Figure [l] and identified a number of 
flux densities from the literature. All flux densities which 
we believe to be reliable are listed in Table JH Flux den- 
sity measurements at low radio frequencies are not easy, 
especially for sources of low surface-brightness; those that 
are available for DA 495 need to be discussed individually. 

The area around DA 495 is covered by the 34 .5 MHz sur- 
vey of iDwarakanath fc Udava Shankaij (jl990[ ). which has 
an angular resolution of 26' x 44' at this declination. An 
apparently unresolved source of flux density 35 Jy is de- 
tected, offset from the position of DA 495 by 12f5 in right 
ascension (to the west) and 5' in declination (to the north). 
The offset may arise from compact sources included in 
the broad beam (see below). The nearby SNR G65. 1-1-0.6 
( Landecker et al.lll99 0) can be identified in the 34.5 MHz 
image. Its expected flux density at 34.5 MHz, extrapo- 
lated from higher frequency measurements, is of the order 
of 50 Jy, but its surface brightness is lower than that of 
DA 495 because of its larger size (1.5° x0.8°). We consider 
the 34.5 MHz data to be generally reliable in this area, and 
adopt an upper limit of 35 Jy for possible emission from 
DA 495 at 34.5 MHz. This is ~7 times the rms sensitivity 
of 5 Jy quoted for the survey. 
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The best low-frequency data co me from the new V LA 
Low-Frequency Sky Survey (VLSS, ICohen et aLll2007t) at 
74 MHz. We obtained an image of the area surrounding 
DA 495 from the VLSS data base@- The angular resolu- 
tion is 80", by far the best of any of our data points below 
327 MHz, and fully adequate to isolate the small-diameter 
sources in the field. The survey is sensitive to structures 
up to ~1° in size, larger than DA 495. The survey noise is 
100 mJy/beam, but in the area around DA 495 the noise is 
higher, about 150 mJy/beam, because of the proximity to 
Cygnus A. DA 495 is not visible in the 74 MHz image, not 
even after smoothing to a resolution of 3'. To obtain an 
upper limit for the flux density at this frequency, we con- 
volved our 1420 MHz image to this resolution, and added 
a scaled version to the 74 MHz image. The scaled image 
was just detectable when its integrated flux density was 
14 Jy, and we adopted this value as the upper limit of 
the 74 MHz flux density of DA 495. We tested this upper 
limit by adding our extrapolated 1420 MHz image to ar- 
eas around the position of DA 495 as well. This convinced 
us that we were not artificially lowering the upper limit 
by adding the 1420 MHz image to 74 MHz emission from 
DA 495 that is just slightly below the detection threshhold. 

Flux densities for DA 495 at 83 and 1 1 1 MH z appear 
in the list published by iKovalenko et al.l (|l994l ). These 
observations were made with fan beams of extent a few 
degrees in the N-S direction. They should perhaps also be 
regarded as upper limits because of the po ssibility of confu- 
sion in the wide fan beams. We note that I Kovalenko et al.l 
(|1994[ ). while listing flux densities for DA 495, give only an 
upper limit for the flux density of G65. 1+0.6. 

Taking the measured flux densities in Table [2] at face 
value we derive the integrated spectrum shown in the left 
panel of Figure [31 fitted reasonably well over the entire 
frequency range by a spectral index of a = —0.59 ± 0.1. 
However, the compact sources in the field around DA 495 
and the superimposed H II region contaminate this spec- 
trum. The compact sources are probably extragalac- 
tic and should m ostly have sp ectral indices steeper than 
a — —0.59 (e.g. ICondonI [1984^ : they will consequently 
contribute strongly to the low-frequency flux densities. 
The H II region will significantly contaminate the high- 
frequency fiux densities because its spectral index is flatter 
than that of the non-thermal emission. We now proceed 
to derive an integrated spectrum for DA 495 corrected for 
these extraneous contributions. 

We have obtained the flux densities of the three most 
prominent compact sources in the field shown in Figure [T] 
at frequencies between 74 MHz and 10550 MHz, using our 
new data and values from the literature. Flux densities 
of these sources are listed in Table [3] and their positions 
and fitted spectra in Table S) The 74 MHz flux density 
for Source 1 is taken from the VLSS point source cata- 
log; the flux densities for Sources 2 and 3 were derived 
by integration from the 74 MHz image. Flux densities for 
DA 495 after correction for compact sources are shown in 
the fourth column of Table [2 

We expect the H II region to affect fiux densities at high 
frequencies, because its spectral index should be a ~ —0.1 
for optically thin thermal emission. The contribution from 
the H II region was easily measured using our images at 



1420 MHz and 10550 MHz. After convolving to an angu- 
lar resolution of 2f45, we plotted brightness temperature at 
1420 MHz against that at 10550 MHz (using images from 
which compact sources had been subtracted) selecting only 
data from the north-east section of the PWN which is un- 
likely to be contaminated by H II region emission. The 
comparison implied a spectral index of —0.85 for this part 
of the PWN. Assuming that this is the spectral index of 
the non-thermal emission between 1420 and 10550 MHz, 
and that the H II region has a spectral index of —0.1, we 
made a point-by-point separation of emission into ther- 
mal and non-thermal contributions. The result is shown 
in Figure [21 The non-thermal emission has strong circu- 
lar symmetry around the central depression, and closely 
resembles our 408 MHz image, where we expect the con- 
tribution of the thermal component to be negligible. The 
thermal component is confined to the south-west, and co- 
incides closely with the infrared source, as expected. Be- 
cause of the similarity of the non-thermal map in Fig. [2l 
and the 408 MHz image in Fig. [1] the coincidence of the 
location of the thermal source in Fig. [21 and the infrared 
source in Fig. [H and the successful separation, we can be 
confident that our assumption of a constant spectral in- 
dex for both the thermal and the non-thermal emission 
components is appropriate. 

The contribution of the H II region to the fiux density 
of DA 495 is 440 ± 100 mJy at 1420 MHz. We used this 
value with a spectral index of —0.1 to calculate the con- 
tribution of emission from the H II region at frequencies 
of 1420 MHz and higher, and subtracted the calculated 
value to obtain an integrated fiux density for the emission 
from the PWN alone. At lower frequencies we did not take 
the emission from the H II region into account, since the 
contribution of the thermal emission should become negli- 
gible and we cannot determine the frequency at which the 
thermal emission becomes optically thick. Flux densities 
for DA 495 after correction for the contribution of the H II 
region are shown in the fifth column of Table [H 

The right panel of Figure [31 shows the integrated spec- 
trum of DA 495 after correction for both compact sources 
and the superimposed H II region. The corrected data 
strongly suggest a break in the spectrum of DA 495 at 
about 1 GHz. Fitting spectra to frequency points below 1 
GHz yields a spectral index of a = —0.45 ± 0.10. Fitting 
to fiux densities above 1 GHz yields a = —0.87 ± 0.10. 
We consider the spectrum above 1 GHz to be well deter- 
mined with a two-parameter fit to 6 well defined fiux den- 
sity measurements. Below the break we have seven mea- 
sured fiux densities two of which are upper limits. There 
are possible errors in the spectral index below the break 
because of the poor angular resolution in that frequency 
range. The quoted error (± 0.10) is the formal Icr-error 
in the fit, but there are probably larger systematic errors 
which we have difficulty estimating; somewhat arbitrarily, 
we adopt a probable error of 0.20. 

These reservations notwithstanding, we believe that this 
analysis has demonstrated convincingly that DA 495 has 
a spectral break at 1.3155:2 GHz of Aa = 0.42 ± 0.22. 
The spectral break may be larger than 0.4, since the low- 
frequency fiux densities are probably upper limits, not 
definite measurements. Further, we have subtracted the 
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contributions of only three compact sources. Some contri- 
bution from such sources will remain in our data, so that 
the plotted values are again to some extent upper limits. 
The location of the break frequency, however, should be 
well defined. Its determination is dominated by the higher 
frequency fluxes of the low frequency part of the spectrum. 
These values are well constrained, while some of the low 
frequency fluxes are upper limits. These affect only the 
slope of the spectrum below the break. 

3.3. Polarization images and rotation measure 

Figure [4] shows our polarization images for the area 
around DA 495, t ogether with data f rom the Effelsberg 
2695 MHz survey (jDuncan et al.lll999r ). At all frequencies 
the polarized emission is confined to a roughly circular re- 
gion about 15' in diameter. No polarized emission is seen 
corresponding to the south-west extension, supporting the 
identification of this as thermal emission. No depolariza- 
tion effects can be associated with this thermal emission 
either, and the H II region which generates it is likely more 
distant than DA 495. 

Integrated over the entire source, the fractional polar- 
ization at the higher frequencies is about 25% (Table [5]). 
Even at 1420 MHz the integrated polarization has dropped 
only to 12%, still half the value at the higher frequencies. 
At 1420 MHz a RM variation of only 36 rad/m-^ is required 
to rotate the polarization angle by 90°, and superposition 
of two signals whose polarization angles differ by 90° pro- 
duces a net signal with no apparent polarization. (At our 
other frequencies RM variations of about 130 rad/m^ at 
2695 MHz, 400 rad/m^ at 4850 MHz, and 1900 rad/m^ at 
10550 MHz are required to generate a Faraday rotation of 
90°.) This implies that there cannot be any rapid changes 
in rotation measure on scales smaller than the 1420 MHz 
beam, because otherwise the 1420 MHz emission would be 
beam depolarized. 

At the two highest frequencies the polarized intensity 
has a remarkable bipolar distribution, seen particularly 
well in the 4850 MHz image (Fig. 2]). Inspection shows 
that the bipolar structure is centered on the pulsar, im- 
mediately suggesting that the polarization structure is de- 
termined by the geometry of the magnetic field, itself tied 
to the pulsar. We see two prominent lobes of polarized 
emission, north-east and south-west of the pulsar. At the 
three high frequencies, from 10550 to 2695 MHz, the lo- 
cation of these lobes seems to rotate systematically clock- 
wise, in the sense opposite to the rotation of polarization 
angle that arises from the positive RM (see below). Sym- 
metrically on either side of the central ridge at 4850 and 
10550 MHz there are small regions of low polarized inten- 
sity. The total intensity remains high at these locations, 
and we interpret the low fractional polarization as an indi- 
cation that polarized emission generated at deeper layers 
within the source is superimposed on emission generated in 
closer layers with a significantly different polarization an- 
gle; vector averaging reduces the net polarization. This is 
largely an emission effect rather than a Faraday rotation 
effect because it is seen at high radio frequencies where 
Faraday rotation is low. 

With linear polarization observations at four frequen- 
cies we should be able to investigate the structure of RM 
across the source. To this end, all observations were con- 



volved to the lowest resolution, 4f3, that of the 2695 MHz 
measurement. However, it proved impossible to calculate 
a RM map: at many points the calculation failed to fit 
a slope to the four measured polarization angles within 
acceptable error limits. The probable cause is "Faraday 
thickness" at the lower frequencies, an effect analogous 
to optical thickness. Depth de polarization - also known 
as differential Faraday rotation (lBurnlll966t ISokoloff et al.l 
il998il - causes the emission from the deeper levels of the 
source to be totally depolarized at low frequencies, so that 
the polarized emission we observe comes only from the 
nearer layers of the source. Faraday thickness seems to af- 
fect the polarization images at 2695 MHz and 1420 MHz. 

Confining our attention to 4850 and 10550 MHz, we cal- 
culated a map of RM at a resolution of 4f3, displayed in 
Figure [51 RM varies from 200 to 250 rad/m^ on the two 
major emission lobes, north-east and south-west of the 
pulsar, to as much as 400 rad/m^ in the structures to the 
east and west. To identify the Faraday thin and Faraday 
thick areas we calculated the expected polarization angles 
at 2695 MHz by extrapolating from the observed angles 
at 4850 and 10550 MHz and compared the results with 
values actually observed at the lower frequency. The po- 
larization angle in the south-western lobe is as predicted, 
while the angle in the north-eastern lobe is about 50° from 
the predicted value. Observed polarization angle in two 
patches lying approximately 4' east and west of the pul- 
sar (these patches are prominent in the 10550 MHz image 
of Figured]) is about 90° from the predicted value. Fara- 
day thickness is clearly significant in this object. We also 
calculated a RM map between 4850 and 10550 MHz at 
the higher resolution of 2 '45 to examine RM fluctuations 
on smaller scales. The RM maps at the lower and higher 
resolution look virtually identical, confirming our impres- 
sion that there are no strong RM fluctuations on small 
scales. We also calculated intrinsic polarization angles at 
the highest resolution and derived a map of the magnetic 
field (projected to the plane of the sky). This is also shown 
in Fig. El 

Even though DA 495 is Faraday thick at frequencies 
of 2695 MHz and below, we calculated a RM map from 
the four bands of the CGPS 1420 MHz polarization data 
(Fig. [5]). Surprisingly, these values of RM, which probe 
only the nearer layers of the source, are of the same or- 
der {^ 200 rad/m^) as the values derived between 4850 
and 10550 MHz where the entire line of sight through the 
PWN is contributing. If RM was positive right through the 
source the values at high frequencies should be higher than 
the values near 1420 MHz. This indicates that there must 
be regions of negative RM in the deeper layers whose ef- 
fects are canceled by positive rotation nearer the front sur- 
face. A possible field configuration that could produce this 
is a dipole field with the dipole axis pointing away from us. 
Almost the entire front part of the nebula would contain 
magnetic field lines that point towards us, explaining why 
we observe only positive rotation measure at 1420 Mhz. 
However, deeper inside the nebula we would find mag- 
netic fleld lines pointing away from us, and the Faraday 
rotation produced by that part of the field would partially 
cancel rotation produced in the near part. The exact RM 
configuration seen would depend on how deep we can look 
inside the nebula at 1420 MHz or other frequencies. 
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4. DISCUSSION 
4.1. Is DA 495 a PWN or a Shell SNR? 

In the introduction we reviewed the arguments for the 
classification of DA 495 as a PWN. In brief, the emissiv- 
ity is centrally concentrated and drops smoothly to zero 
at the outer edge without any sign of a steep shell, indi- 
cating a central energy source. This c onclusion is corrob- 
orated by X-ray observations (jArzoumanian et al. 2004. 
l2008t ) that show a cei itral compact object su rrounded by a 
non-thermal nebula. IVelusamv et al.1 (|l989t ) suggest that 
DA 495 could be a shell remnant with an unusually thick 
shell. While this is a conceivable interpretation, there is no 
evidence of a steep outer edge to the shell and the central 
compact source discovered in X-ray is not located inside 
the radio hole, but about 2 ' away from it suggesting that 
the radio hole is not in the centre of the nebula. 

This paper provides further evidence. The structure in 
polarized intensity is bipolar with mostly radial B- vectors 
at the outer edge, while for a shell-type remnant we would 
expect a shell structure in polarized emission with tangen- 
tial B-vectors. A young shell-type SNR would have ra- 
dial B-vectors, but integrated polarization would be very 
low and a young remnant would have bright X-ray emis- 
sion coming from the expanding shell. The break in the 
integrated spectrum of DA 495 also provides strong evi- 
dence that the object is a PWN. Every known example 
of such a nebula exhibits such a break, while virtually 
all shell remnants have spectra characterized by a single 
power law. While spectral breaks have been claimed for 
some shell remnants, these claims are gradually melting 
away in the face of new data with improved sensitivity 
to extended emission. The only SNR where such a break 
has been convincingly established is S 147 (G180. 0—1.7) 
which has a break in i t s integrated spect rum at ^1.5 GHz 
(|Fiirst fc Reich Ill98d IXiao et al.ll2008D . S147isavery 
old SNR, probably well into the isothermal stage of its 
evolution; it bears no resemblance to DA 495. 

The only tenable conclusion is that DA 495 is a pulsar 
wind nebula, and we will proceed on that basis. 

4.2. Distance, Size, and Emission Structure 

A d istance to DA 495 was measured by iKothes et al.l 
(|2004[ ) on the basis of absorption by H I of the po- 
larized emission from the PWN. The distance obtained, 
d ~ 1.5 kpc, was derived kinematically using a flat rota- 
tion model for the Milky Way with a Galactocentric radius 
of Rq = 8.5 kpc for the Sun and its velocity vq = 220 km/s 
around the Galactic center. The latest measurements 
of the Sun's galactocentric distance give Rq « 7.6 kpc 
which is used in the new dista nce determination method of 
[Foster fc MacWilliamsl (|2006l . and reference therein) . This 
method uses a model for the spatial density and velocity 
field traced by the distribution of H I in the disk of the 
Galaxy to derive a distance-velocity relation that can be 
used to calculate a distance to any object with a known 
systemic velocity. The H I absorption measurements in- 
dicate for DA 495 a systemic velocity of about H-12km/s. 
Two distances, about 1 and 5 kpc, correspond to this ve- 
locity, one closer than the tan g ent p oint and the other 
beyond it. Since iKothes et al.l (|2004f ) found no absorp- 
tion of the PWN at the tangent point and DA 495 ex- 



hibits low foreground absorption in the X-ray observations 
(jArzoumanian et all |2004 I2008D we adopt the closer dis- 
tance as the more likely result. We proceed to use a dis- 
tance of 1.0 ± 0.4kpc to determine physical characteristics 
of the PWN. The error contains the uncertainty in the fit 
and a velocity dispersion of 4kms~^. 

To determine a reliable value for the radius we fit- 
ted a three-dimensional model of emissivity to the total- 
intensity obser vations at 1420 MHz. Our me thod is similar 
to that use d bylLandecker fc Caswell (1X9831), a method de- 
veloped bvlliil (fT967h . While iLandecker &: Caswell (|T98l 
centered their calculation on the depression near the mid- 
dle of the nebula, we centered our ne w calculation on the 
positio n of the pulsar, as measured bv I Arzoumanian et al.l 
). In the center there is an emissivity depression with 



a radius of about 2', and so a physical radius of 0.6 pc. 
Outside, the emissivity is nearly constant to a radius of 7'.5, 
where it drops abruptly, al most to zero. This translates 
to an outer radius of 2.2 pc. ILandecker fc Caswelll (|1983( ) 
found that the emissivity increased from the outer edge to- 
wards the centre. This difference can be explained by the 
different location assumed for the centre and the contam- 
ination by the H II reg ion (removed from ou r calcu lation 
b ut included in that of Landecker fc Caswelll (|1983D ). 

ILandecker fc Caswelll (|l983l) assumed that the promi- 
nent depression marked the center of the PWN, but we 
know now that the pulsar lies at one edge of that depres- 
sion (jArzoumanian et al.l |2004 I2008D . not at its center. 
One might suppose that the pulsar has moved since it was 
formed, but we see from Figs. |3] and [5] that the structure 
in polarized intensity, and so the magnetic field structure, 
is centered on the pulsar. Since the radio appearance of 
the nebula is dominated by particles injected early in the 
life of the pulsar (see Section 4.3), it appears that the pul- 
sar has not moved significantly, and we must find another 
explanation. 

In order to achieve higher angular resolution at 
1420 MHz, we merged the 1420 MHz map from Fig.[T] 
with data fro m the NRAO VLA Sky Survey (NVSS; 
ICondon et anil998) . The NVSS has a resolution of about 
45". The resulting map is shown in Fig.jS) Inspec- 
tion of Fig.jS] shows evidence of a second less obvious de- 
pression in total intensity, at RA(2000) = 19'' 52™ 33^ 
DEC(2000) = 29° 26'. This eastern depression is further 
from the pulsar than the more obvious, and deeper, hole 
to the west, and the nebula appears to extend to slightly 
larger radius beyond the eastern hole. Nevertheless, we 
regard the eastern depression as a significant feature. In 
order to produce these two emissivity depressions, there 
must be a region around the pulsar which is relatively free 
of magnetic field and/or particles. This brings to mind the 
Crab Nebula's "synchrotron bays" . iFesen et al.l (119921 ) as- 
cribe them to an equatorial torus of material ejected by 
the progenitor star before it exploded. They suggest that 
the magnetic field of the nebula is wrapped around this 
torus preventing relativistic particles from entering it. A 
saddle-shaped depression in total power emission that con- 
nects the two holes in DA 495 (see Fig. [5]) supports a torus 
as an explanation for these features. In the Crab Neb- 
ula the synchrotron bays are located at the edge of the 
emission region, but DA 495 is probably much older and 
the nebula has had time to flow around the torus. The 
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distance between the two holes is about 5' translating to 
1.45 pc at a distance of 1 kpc. This indicates an average 
expansion velocity of about 70 ^0"°" y^ km/s, where t is 
the age of DA 495. For an age of 10000 yr this expansion 
velocity would be about 10 % of that of the Crab Neb- 
ula's synchrotron bays, which are expanding at about 900 
and 600 km/s for t he western and eastern bay respectively 
(JFesen et al.l 119921 ). We cannot completely rule out that 
movement of the pulsar is responsible for the off-center 
location of the pulsar between the two holes, but a differ- 
ence in expansion velocity between the two sides, similar 
to that observed in the Crab Nebula, could explain this 
asymmetry. 

4.3. The Nature of the Spectral Break 

To investigate the nature of the break in the radio spec- 
trum of DA 495 we first estimate the nebula's age by com- 
paring its energetics with the properties of other PWNe 
of known age. We assume for the moment that the pul- 
sar has a symmetrical pulsar wind and that equipartition 
between particle and magnetic field energies has been es- 
tablished. We follow the analysis made by iKothes et al.l 
((2006a) of the PWN G106.65-h2.96, although, unlike that 
case, we have no information on the characteristic age of 
the DA 495 pulsar, which has not been detected. Nev- 
ertheless, we can get an esti mate of the current energ y 
loss rate of the pulsar from lArzoumanian et all (J2008D : 
-Bda495 ~ 10^^ erg/s. We now confine our attention to 
pulsars in historical SNRs for which we know the exact 
age, the Crab Nebula (SN 1054), 3C 58 (SN 1181), and 
Gil. 2-0.3 (SN 386) (lGreenl l2004D. We al so add three 
pulsar wind nebulae for which 'K othes et all (|2006aF ) have 
calculated the real ages from the location of the cooling 
breaks in their radio spectra. We extrapolate their evo- 
lution to ii'DA495 to determine possible properties of the 
pulsar in DA 495. 

The evolutio n of the energy loss rat e, E, of a pulsar is 
defined by (e.g. lPacini fc Salvati|[T973l ): 
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where (3 is related to the nature of the braking torque 
(/3 = ^r^T ' ^"^^ "■ ^^ ^^^ braking index of the pulsar) , t is 
the real age of the pulsar, Eq is its energy loss rate at the 
time it was born, and tq is its intrinsic characteristic age 
(t'o = -^ , where Pq is the period of the pulsar at the time 
it was born). Using the information about the current 
characteristic age r and E f or the six pulsars taken from 
the ATNF Pulsar Catalogue (iManchester et al.ll2005H1 and 
their real age we can derive their initial parameters, tq and 
Eq . For the Crab Pulsar and the pulsar i n Vela X the brak - 
ing indic es are know n, yielding a /? of 2.3 (jLvne et al.lll988[ ) 
and 6.0 (jLvne et al.l l 19961. respectively. For the other four 
pulsars we assume a pure dipole field, implying a /3 of 2. 
For Vela X the age depends on the assumed value of /3, 
since we calculate it from the spectral break frequency. 
The resulting intrinsic pulsar properties are listed in Ta- 
ble [71 To indicate how much our calculations depend on 
this factor we used different values of f3 as indicated in 
Table [71 Using equation 1 and taking the energy loss rate 
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of the pulsar in DA 495 as iJDA495 ~ 10'^^ erg/s, we can 
calculate the current age, iDA495, of the nebula and the 
current characteristic age, tda495, of its pulsar using the 
initial values for each of the six pulsars. The results are 
given in Table [51 For the Crab and the Vela X pulsars 
we extrapolated the energy loss rate with two different /3 
values as indicated in Table [51 Whatever basis we adopt 
for the calculation, the age of DA 495 exceeds 20,000 yr, 
which makes this nebula an old object. This derived age 
is valid regardless of the mechanism responsible for the 
break in the radio spectrum. All known PWNe detected 
at radio wavelengths contain energetic pulsars, so the low 
energy loss rate of the DA 495 pulsar is in all likelihood 
the result of age, not of intrinsically low energy. 

We will now proceed by examining the different possi- 
bilites for a spectral break in the synchrotron spectrum of a 
pulsar wind nebula. Virtually all PWNe exhibit breaks in 
their spectra, at frequencies ranging from X-rays to a few 
GHz. Table [6l lists the break frequencies that have been 
identified to date. The spectral break in DA 495 occurs 
at 1.3 GHz, significantly lower than for any other pulsar 
wind nebula. A spectral break in a PWN can arise in 
two ways: it may reflect a break in the injected electron 
spectrum, and so be a mostly pulsar-dependent effect, or 
it may arise from synchrotron cooling, in which case the 
break frequency depends on the magnetic field strength 
and the age of the relativistic electrons in the nebula. Of 
course, both mechanisms should be at work in all PWNe. 
The suspected nature of each spectral break is indicated 
in Table [51 

We first consider the possibility that the break in the 
radio spectrum results from a break in the injected elec- 
tron spectrum. Such breaks in the spectra of PWNe are 
known (e.g. in the Crab Nebula and 3C 58), but they are 
imperfectly understood. Since the electron population in 
the nebula is dominated by early injection, the break must 
have been created early in the life of the pulsar. The en- 
ergy at which the break occurs depends on conditions in 
the area in which electrons are accelerated. Its location in 
the electron spectrum remains relatively fixed throughout 
the life of the nebula. In contrast, a break due to cooling, 
starts at hig h energy and m oves with time towards lower 
energies (e.g. IWeiler fc Pan agia 1980) ultimately perhaps 
passing by the break produced by the injected spectrum. 
If the magnetic field decays with time, through expansion 
or other causes, both kinds of break will migrate towards 
lower frequencies in the radio spectrum. Their relative 
positions will not be affected. 

A likely location for an injected break is at the energy 
where the electron acceleration timescale equals the syn- 
chrotron loss time, which is in essence a kind of cooling 
break. The difference between this situation and the nor- 
mal synchrotron cooling break is that the number of high 
energy electrons is never high at any stage. Thus, the 
synchrotron spectrum should steepen by much more than 
the 0.5 typical of normal synchrotron cooling, which is ac- 
tually the case for the intrinsic breaks identified in other 
PWNe. For DA 495 we found that the spectrum steepens 
by about 0.4. 

There are two other arguments that make an in- 
jected break an unlikely interpretation for the break in 
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the radio spectrum of DA 495. First, the X-ray syn- 
chrotron spectrum has a spectral index of —0.63 ± 0.27 
([Arzoumanian et al.ll2008l calculated from the photon in- 
dex T: a = — (r — 1)), flatter than the synchrotron radio 
spectrum above the break (—0.87). This is impossible if 
the radio break arises from an injected break. The only 
mechanism that can affect the synchrotron spectrum be- 
yond an intrinsic break is cooling or some other energy 
loss effect, which can only steepen the spectrum but can- 
not flatten it. A new round of particle acceleration late in 
the life of the PWN, an unlikely event, would affect low- 
energy and high-energy electrons alike. Second, there is 
no injected break frequency known that is even remotely 
as low as the 1.3 GHz we have found in DA 495. 

We now consider the possibility that the break at 
1.3 GHz is caused by synchrotron cooling. We will assume 
that the magnetic field throughout DA 495 is relatively 
constant or at least that there is a dominating magnetic 
field creating the cooling break. There are several lines of 
evidence supporting this assumption. First, our calcula- 
tion of the synchrotron emissivity from the total-intensity 
images shows that emissivity is nearly constant in the ma- 
jority of the nebula (Section 4.2). Second, the break in 
the total spectrum is very sharp and there is no signifi- 
cant variation of spectral index across the source (DA 495 
has the same appearance at all radio frequencies, see also 
section 3.2: the separation into thermal and non-thermal 
components). Whatever the cause of the spectral break, 
it is the magnetic field that translates a break in the elec- 
tron energy spectrum to a break in the emission spectrum. 
If the magnetic field varied significantly, the break would 
be smeared across a wide range of frequency; the break 
frequency has a very strong dependence on field strength 
(see Equatio n (2) below). 

IChcvaliei] (|2000D found that the cooling frequency, Vc, 
corresponding to 7 (the Lorentz factor of the electron) at 
which the electrons are able to radiate their energy over 
the age of a PWN, is defined by 

i/c[GHz] = 1.187 • B-^[G] ■ t-^[yr], (2) 

where B is the magnetic field strength inside the syn- 
chrotron nebula. Using the break frequency of Vc ~ 
1.3 GHz and the various interpolated ages, tDA495, we can 
calculate i?req, the required magnetic field inside the neb- 
ula (Table [5]). Integrating equation (1) over the age of 
the nebula gives the total energy inside it and, assuming 
equipartition and a circular structure, we can calculate the 
magnetic field, Bmax- We denote this field value iJmax be- 
cause we have not accounted for any energy loss in the 
course of the evolution of the PWN. Comparing the re- 
sults based on the six pulsars listed in Table [SJ we see that 
it is only in the calculation based on the Crab pulsar and 
the pulsar in Vela X that the maximum magnetic field ex- 
ceeds the magnetic field required for the spectral break. If 
our assumption is valid, then the pulsar in DA 495 must 
have been very energetic with a high energy loss rate at 
the time of formation, like the Crab pulsar or the pulsar 
in Vela X. Given the age of the Vela supernova remnant 
and the fact that i?max is significantly higher than Sreq 
for any /3 it is more reasonable to assume that the pulsar 
in DA 495 is an older version of the Vela pulsar. If we 
shrank the Vela X nebula with its current energy content 
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to the size of DA 495, its cooling break at 100 GHz would 
decrease to about 2.5 GHz due to the resulting increase in 
the magnetic field strength. A few thousand years later 
the break frequency would reach 1.3 GHz. 

Even though we cannot completely rule out the possi- 
bility that the break in the radio spectrum of DA 495 at 
1.3 GHz is caused by a break in the injected electron spec- 
trum, all the available evidence indicates that synchrotron 
cooling is the most likely explanation, and that DA 495 is 
probably an aging Vela X nebula. In this case DA 495 is 
about 20,000 yr old and contains a strong magnetic field 
of about 1.3 mG. This can easily account for the high ro- 
tation measure we have observed. 

4.4. The magnetic field structure inside DA ^95 

The structure of the PWN in total intensity is very reg- 
ular, especially once the confusing H II region is removed 
(see Figure[2|). However, in polarized intensity the appear- 
ance is very different, and at the high frequencies DA 495 is 
very structured. Percentage polarization, even at the high 
frequencies, varies over a wide range across the source. 
The challenge faced in this section is to deduce something 
about the structure of the magnetic field inside DA 495 
and to explain as many of the observed characteristics as 
possible. While we cannot uniquely determine the field 
configuration we show that one relatively simple and plau- 
sible model, a dipole field with a superimposed toroidal 
component, is consistent with the observations. 

4.4.1. The foreground Faraday rotation 

In Figure \7\ we have plotted the foreground RM ob- 
served to all pulsars within 15° of DA 495 as a func- 
tion of distance from us in kpc. The pulsar data were 
taken from the ATNF Pulsar Catalog (|Manchester et all 
[2005 rl. All pulsars within about 5 kpc have negative RM 
values, except for one for which RM = +7 rad/m ; we 
will disregard this exception because it lies at a Galac- 
tic latitude of -1-14°. The negative RM values are con- 
sistent with the known local magnetic field orientation, 
away from the Sun. Pulsars more distant than 5 kpc 
have positive RMs, as expected from the well-known field 
reve rsal between the local and Sagittarius spiral arms 
(e^ Simard-Normandin fc Kronbergjll979l : iLvne fc SmithI 
19891). Averaged over the nearby pulsars the line-of-sight 
component of the magnetic field is around —2.5 fj,G. Using 
the m odels oftTaylor & Cordes (1993) and lCordes fc Lazig 
(|2002l ) for the distribution of free electrons in the Galaxy, 
and the distance to DA 495 of 1.0 kpc we determine 
the average foreground dispersion measure to be DM w 
15 cm~'^ pc. With an average foreground magnetic field of 
—2.5 fiG parallel to the line of sight this produces a fore- 
ground RM of —30 rad/m^. The positive RM values deter- 
mined from the emission of DA 495 must arise from Fara- 
day rotation within the PWN itself or occur in a discrete 
foreground magneto-ionic object which does not itself pro- 
duce enough emission to be detected in our observations. 
Since we find that almost all of the polarized emission from 
DA 495 is Faraday thick at 2695 MHz, most of the positive 
RM must originate in DA 495, because Faraday thickness 
can only be produced if the rotation is generated in the 
same volume elements as the polarized emission itself. A 
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foreground object can merely rotate or beam depolarize 
the entire emission coming from behind it. 

4.4.2. Magnetic field configurations inside DA 495 

The magnetic field distribution seen in Figure [5] is 
the projection onto the plane of the sky of the three- 
dimensional field within DA 495, as deduced from our ob- 
servations at 4850 and 10550 MHz on the assumption 
that the source is Faraday thin at these frequencies. The 
first impression is of a highly regular elongated magnetic 
field configuration with field lines following a central bar 
and spreading out radially to the north and south. This 
strongly resembles the inner part of a dipole field projected 
onto the plane of the sky. The field lines must close, pre- 
sumably at large radius. However, in Figure [5] we are see- 
ing only the component perpendicular to the line of sight. 
The relatively low percentage polarization even at very 
high frequencies of about 25 % - intrinsically it should be 
about 70 % - cannot be explained by depolarization caused 
by Faraday rotation, because at high frequencies this ef- 
fect is negligible: a RM of 400 rad/m^, the highest we have 
detected, rotates the polarization angle by a mere 20° at 
10550 MHz. Hence, the depolarization must be caused by 
the superposition of polarized emission generated at differ- 
ent points along the line of sight with intrinsically different 
polarization angles, averaging to a low net polarized frac- 
tion. In a pure dipole field the component perpendicular 
to the line of sight, which determines the observed polar- 
ization angle, does not change significantly along different 
lines of sight through the nebula. There must be an addi- 
tional magnetic field present that also creates synchrotron 
emission but has an intrinsic magnetic field configuration 
perpendicular to, or at least very different from, the dipole 
field. A toroidal field component wrapped around the axis 
of the dipole would match that description. Toroidal field 
struc tures have been observed in Vela X (iDodson et al 
200l and in the PWN in SNR G106.3-H2.7 (jKothes et a" 
2006a[) . although most other pulsar wind nebulae show 



only an elongated radial field, indicating the presence of a 
dipole. The patches of polarized emission to the east and 
west of the pulsar have a magnetic field structure perpen- 
dicular to that expected for a pure dipole field and could 
be a hint of a toroidal magnetic field. The two promi- 
nent "holes" in polarized intensity would then mark areas 
where the polarized emission generated by the dipole com- 
ponent and that generated by the toroidal component are 
equally strong but have orthogonal polarization angles so 
that the net emission is intense but unpolarized. Even 
though both magnetic field configurations may be present 
it is clear that the dipole field is the dominant one for emis- 
sion, as seen in the magnetic field projected to the plane 
of the sky, and for Faraday rotation. If most of the RM 
was produced by the toroidal field one half of the nebula 
would show a negative RM and the other a positive one. 

One possible magnetic field configuration that could 
explain the radio polarization observations of DA 495 is 
shown in Fig. [51 This is no more than a naive sketch, but 
it fits the available observations remarkably well. A sim- 
ilar combination of dipole and toroidal field components 
can explain the polarization appearance of other PWNe; 
the polarized emission structure, the magnetic field vectors 
projected to the plane of the sky, the fractional polariza- 



tion, and the Faraday rotation structure all depend on the 
viewing a ngle and which one of the two components is 
dominant ([Kothes et al.ll2006af ). 

4.4.3. The apparent rotation of the polarization structure 

One curious feature of the radio emission from DA 495 
is the seemingly systematic rotation with frequency of al- 
most the entire polarized emission structure. At the three 
higher frequencies we can see the two prominent lobes ro- 
tating from north-east and south-west at 10550 MHz to 
west and east at 2695 MHz. Near the center the two po- 
larization holes (or, equivalently, the emission ridge) also 
rotate clockwise. At 2695 MHz we can no longer see this 
rotation near the center, probably because of the larger 
beam. This effect must be the product of depth depolar- 
ization in a highly ordered magnetic field structure. 

A possible explanation is the following. At the intrin- 
sic (zero wavelength) location of the polarization holes the 
magnetic field lines of the toroidal and dipole fields are 
perpendicular, and the polarization averages out. Close 
to the center of the PWN the dipole field lines are almost 
all parallel, while the toroidal field lines change with az- 
imuthal angle within the torus. If the Faraday rotation of 
emission generated in the toroidal field differs from that 
generated in the dipole field then, at lower frequencies, this 
averaging no longer occurs at the "intrinsic" location but 
at a location where the dipole and toroidal field lines are at 
some angle other than 90°. The location of the holes then 
moves systematically with frequency, in the sense opposite 
to Faraday rotation, as observed. 

If the apparent rotation of the polarized structure is a 
Faraday rotation effect, then it should show a A^ depen- 
dence, and we can test this. We define ip as the position 
angle of a line through the centers of the lobes {ip increas- 
ing from north to east). For 10550 MHz i/j — 45° and for 
4850 MHz ip = 15°. Drawing a second line through the 
polarization holes, and assuming it is perpendicular to the 
fine through the lobes, we find ip = 38° at 10550 MHz 
and ip = 10° at 4850 MHz. Averaging these estimates, we 
derive a "structural rotation" of about —170 rad/m^ be- 
tween 10550 and 4850 MHz and an intrinsic ip of 50°. As a 
test, we extrapolated to 2695 MHz, and predict ip — —70°, 
which is very close to the —90° we observe. A second 
test comes from comparison with our RM map (Figure 
5). Faraday rotation seems to be dominated by the dipole 
field, and the RM should then show a maximum or mini- 
mum at the center of the lobes because in those directions 
the angle between the line of sight and the field lines is 
smallest. A line from the RM peak in the north-eastern 
lobe to the pulsar has an angle of 55° with north, compa- 
rable with the intrinsic ip of 50° calculated above. 

Since the spin axis of the pulsar is the only axis of sym- 
metry in the rotating pulsar frame, we can assume that 
it is parallel to the axis of the nebula's dipole field, which 
then would reflect the rotation-averaged pulsar field. We 
conclude that the spin axis of the pulsar, the main axis of 
the toroidal field, and the axis of the dipole magnetic field 
component all have an angle ip « 50° projected onto the 
plane of the sky. Further, the viewing angle, 0, between 
the line of sight and the equatorial plane of the dipole (see 
Fig. [ID, must be larger than 0°. At 6* = 0° the toroidal field 
component would be constant projected onto the plane of 
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the sky and we would not observe the structural rotation 
effect. 

4.4.4. Orientation of the dipole field 

In this section we develop a simple model of the mag- 
netic field configuration in DA 495. We test the model by 
computing RM along a few selected lines of sight through 
the nebula, and achieve plausible agreement with measure- 
ments. We have not attempted detailed modeling of the 
field. 

At a distance of 1 .0 kpc the diameter of DA 495 is 4.4 pc; 
we assume this is the maximum path length through the 
nebula. We deduced an internal field of about 1.3 mG from 
our comparisons with the Vela X nebula (Section 4.1 and 
Table [5]) and we assume that the field is constant at this 
value throughout the nebula. We know the projected field 
distribution (Figure [5]) and we know that the dipole field 
should be the dominating factor for the internal Faraday 
rotation. For these calculations we assume a dipole field 
as displayed in Figure [H We confine our attention to the 
main axis of the dipole, because any toroidal component 
will be perpendicular to the line of sight and will not con- 
tribute to RM (and so we ignore the complication of the 
toroidal component). Magnetic field strength and syn- 
chrotron emissivity are assumed constant throughout the 
nebula (the evidence for these assumptions is presented in 
Sections 4.2 and 4.3), and the intrinsic polarization angle 
is set to be constant at all points along every line of sight 
through the nebula. We calculate the electron density rig 
within the nebula needed to generate the observed RM af- 
ter correction for the foreground. We make this calculation 
for various viewing angles, 9. In our calculations we are 
looking from the right hand side of Fig. [5] into the nebula. 
We consider we have succeeded if we can find an electron 
density within the nebula, assumed constant, and a view- 
ing angle which together provide RM values reasonably 
consistent with the observed RM along different lines of 
sight. The foreground RM of approximately —30 rad/m^ 
is taken as a measure of the acceptable uncertainty. 

Consider first the appearance at 9 = 0°. In the north- 
ern half the emission generated in the deeper levels first 
passes through a magnetic field that is directed away from 
us on the far side of the axis and then through a region 
where it is pointing towards us on the near side. Since 
the emission is generated uniformly along the line of sight, 
the bulk of the emission passes through the magnetic field 
that is pointing towards us and the resulting rotation mea- 
sure is positive. However, in the bottom half the structure 
is reversed, and this should lead to a negative RM. We 
would expect to see a symmetrical RM map, positive to 
the north and negative to the south with peaks of RM of 
equal absolute value. Now consider the appearance at an 
angle 9 > 0° as portrayed in Figure [S] For lines of sight 
above the center of the nebula, positive RM dominates, 
as above. For lines of sight passing below the center the 
emission from the far side of the nebula emerges with a 
strong positive RM because the path is long; there is some 
path through the near side that contributes negative rota- 
tion, and there is some reduction of net RM, but the path 
through the near side is relatively short, so that the final 
RM value remains strongly positive, indeed more positive 
than in the upper lobe. Our observations do show posi- 



tive RM throughout, but with a higher RM value on the 
north-east lobe. Therefore to represent the observed situ- 
ation, the field distribution in Figure[H]must be rotated by 
180° around the line of sight, so that the dipole is directed 
downwards and away from us. 

In order to test whether this simple model can pre- 
dict the observed properties of DA 495, we calculated the 
RM values that would be observed between 10550 and 
4850 MHz for three different paths through the nebula, 
through the center and 1.2 pc above and below the cen- 
ter. At the center we observe a rotation measure of about 
170 rad/m^, which amounts to 200 rad/m^ after correction 
for the foreground RM. At the RM peak on the north- 
eastern lobe, 1.2 pc above the center, we observe about 
310 rad/m^ and at the same distance on the south-western 
lobe we find about 250 rad/m^. We calculated the electron 
density, Ug, required to produce these RM values for the 
three lines of sight through the nebula for varying view- 
ing angles 9. The results are displayed in Fig. [TOl We 
find that the model, although approximate, can reproduce 
the observed RM values with reasonable values for the 
parameters. The small rectangles identify the region of 
satisfactory agreement. Within the limitations of our as- 
sumptions, the electron density in the nebula lies between 
0.3 and 0.4 cm~'^ and the viewing angle 9 is between 33° 
and 40°. 

In our model we assumed constant magnetic field 
throughout the nebula. Even though the total emissiv- 
ity and hence the total magnetic field strength inside the 
nebula must be relatively constant, we know that the mag- 
netic field strength in a dipole field decreases with distance 
from the centre. We use another simple model, going to 
the other extreme, to assess the effect of this assumption: 
we collapse all the rotation measure along the line of sight 
onto the central axis of the dipole. All the emission gen- 
erated beyond the axis is rotated and then added to the 
emission from the near side. Due to the symmetry of the 
dipole, the collapsed RMs at the two positions above and 
below the centre used in our calculations would be the 
same. However, by changing the viewing angle 9 we would 
change the pathlength through the nebula in front of and 
beyond the central axis. For constant emissivity we de- 
rive an angle 9 of 10°, in contrast to the ~36° derived 
above. Under both assumptions, and also under condi- 
tions of varying emissivity within the nebula, the dipole 
must be directed away from us towards the south. 

Finally, in Fig. [TT] we illustrate the three dimensional 
orientation of the dipole field inside DA 495 as determined 
in this section. On the left hand side we display the mag- 
netic field projected to the plane of the sky. The magnetic 
field axis has an angle of about 50° east of north and the 
magnetic field is pointing south-west. On the right hand 
side the dipole field orientation along the line of sight is 
shown, indicating that the dipole is pointing away from us 
towards the south-west. 

4.5. Limits on a surrounding shell-type SNR 

Is DA 495 a pure pulsar wind nebula, like the Crab Neb- 
ula, or does it have a surrounding SNR shell? With our 
new data we were able to search for such a shell. None 
was found, and we can place sensitive upper limits on 
any such emission. The noise in the area around DA 495 
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in our 1420 MHz image, which is the most sensitive we 
have for such a search, is about 40 mK rms. A Sa sig- 
nal would be about 1.0 mJy/beani, and we would need 
this peak flux density for a positive detection. If the in- 
ner edge of this shell is at the outer edge of the PWN 
and its shell thickness is 10%, typical for an adiabatically 
expanding SNR shell, the upper limit for the total flux 
density is about 90 mJy at 1420 MHz. An older rem- 
nant would have higher compression and a thinner shell, 
resulting in an even lower flux density limit. The 1 GHz 
flux density limit is about 100 mJy, using a spectral in- 
dex of a = —0.5, which is a typical value for a shell- 
type SNR. This translates to a maximum radio surface 
brightness of S'i^^^hz = 5-4 x IQ-^^ Watt m'^ Hz"! sr"!, 
slightly lower than the lowest radio surface-brightness 
known, that of the SNR G156.2- I-5.7 which has Ei ghz = 
6 X 10^23 -y^g^tt j^-2 JJ2-1 sr-i (|Reich et al. 1992). 

The 60^m image of the DA 495 field (Fig. [T|) shows in- 
frared emission to the north-east of the SNR. We found H I 
in the CGPS data, which is likely related to this emission, 
because it exhibits a similar shape. This H I, however, 
is found at velocities near the tangent point, far removed 
from the systemic velocity of DA 495. We do not believe 
that any of this infrared emission is related to the SNR. 

5. CONCLUSIONS 

We have imaged the radio emission from the pulsar wind 
nebula DA 495 in total intensity and polarized emission 
across a range of frequencies from 408 to 10550 MHz. Af- 
ter correcting for the flux contribution of extraneous su- 
perimposed sources we have demonstrated that there is a 
break in the spectrum at about 1.3 GHz which we can con- 
vincingly attribute to synchrotron cooling. The spectral 
break at such a low frequency, the lowest known for any 



PWN, can only be produced if the pulsar originally was 
very energetic, with a high energy loss rate. We consider 
that DA 495 is an aging pulsar wind nebula which has a 
pulsar that in its earlier stages resembled the Crab pul- 
sar or the pulsar in Vela X. The magnetic field inside the 
nebula is remarkably well-organized and has the topology 
of a dipole field centered on the pulsar; the field probably 
also has a toroidal component, wrapped around the dipole 
axis. The rotation measure is high, as would be expected 
in a nebula where the magnetic field has the high value 
of '^1.3 mG, and the nebula becomes Faraday thick at a 
very high frequency. We have set a very low upper limit 
for the surface brightness for any shell of SNR emission 
surrounding DA 495. The absence of a shell implies that 
the SN explosion occurred in a low-density environment, 
and this also explains how the magnetic field structure has 
remained so regular to an age of '^20,000 yr. 
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Table 1 
Observation characteristics 



Frequency [MHz] 



408 



1420 



4850 



10550 



Telescope 
Observation Date 



Half-power beam (') 

rms noise, /, [mJy/beam] 
rms noise, PI, [mJy/beam] 
Calibrators, /, (Flux Density, Jy) 



Calibrators Q, U 
Linear Pol. [%] 
Pol. Angle [°] 
Coverages'^ 



DRAG 


DRAG 


Effelsberg 


Effelsberg 


January 


January 


February 


January 
to July 


2002 


2002 


1996 


1990 


2.9 X 5.95 


0.82 X 1.75 


2.45 


2.45" 


(R.A. X Dec.) 


(R.A. X Dec.) 






7 


0.4 


0.5 


<1.5'' 




0.4 


0.3 


0.4 


3C48 (38.9) 


3C48 (15.7) 


3C286 (7.5) 


3C286 (4.5) 


3C147 (48.0) 


3C147 (22.0) 






3C295 (54.0) 


3C295 (22.1) 








3C286 


3C286 


3C286 




9.25 


11.3 


11.7 




33 


33 


33 






14 


5 



" Original resolution 1.'15. Smoothed to 2f45 for presentation. 
** Estimated. There is no region in the image entirely free of emission. 



Effelsberg observations only. 
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Table 2 
Integrated flux densities of DA 495 



Frequency [MHz] 


Beam 


^meas 


[Jy] 


Scs-r [Jy] 


Shii-t 


[Jy] 


reference 


34.5 


26' X 44' 


< 35.0 


<24.8 






1 


74 


1.'33 X 1.'33 






<14 






5 


83 


0.2° X 5.2° 


18 ± 


5 


13.0 ± 5 






2 


111 


0.15° X 3.9° 


17 ± 


5 


13.0 ± 5 






2 


318 


14' X 16' 


9.7 ± 


2.2 


8.1 ± 2.2 






3 


327 


l.'l X 0f6 






7.8 ± 1.0 






4 


408 


2.'8 X 57 






6.5 ± 0.6 






5 


430 


8.'4 X 9f6 


8.7 ± 


4.9 


7.5 ± 4.9 






3 


610.5 


16' X 18' 


7.0 ± 


1.5 


6.0 ± 1.5 






6, 7 


1420 


0.'8 X lf6 






4.0 ± 0.2 


3.6 ± 


0.3 


5 


1665 


22' 


4.4 ± 


0.5 


4.0 ± 0.5 


3.6 ± 


0.5 


7 


2695 


4^3 






2.7 ± 0.2 


2.3 ± 


0.3 


8 


2735 


4.'9 X 57 


2.8 ± 


0.4 


2.7 ± 0.4 


2.3 ± 


0.4 


7 


4850 


2'Ab 






1.59 ± 0.1 


1.2 ± 


0.1 


5 


10550 


lfl5 






1.12 ± 0.1 


0.76 ± 


0.18 


5 



Notes: Smeas is the measured integrated flux density of DA 495. Scs-r is the measured flux density minus the total flux 

density of the three brightest compact sources (see text) appropriately scaled with frequency. For the new data 

presented in this paper, and for other high-resolution measurements, this is the integrated flux density of DA 495 

without the compact sources. Suii-r is the integrated flux density minus an allowance for the H II region which 

ove r laps t he SNR (see text) 



Refer ences: iV) iDwarakanath fc Udava Shankarj (|1990D. (2) Koval cnko etal] (Il994l). (3\ lDickel fc DeNovei) (Il975l) . (4) 



IVelusamv et al.l (Il989f ). fS) present work. r6l lDickelet all (I1971D . m IWillid (Il973f ). (8') lReichet all (ll990D 
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Table 3 
Flux densities of compact sources [mJy] 



Frequency, MHz Source 1 Source 2 Source 3 Reference 

74 1470 ± 180 1200 ± 400 900 ±300 3 & 10 

330 613 ± 50 711 ± 70 328 ± 30 1 

365 557 ± 67 739 ± 52 294 ±28 2 

408 650 ± 80 504 ± 50 200 ±40 3 

1420 139 ± 5 280 ± 8 93 ± 5 4 

1420 131 ±5 244 ± 25 80 ± 6 5 

1420 134 ±8 288 ± 15 83 ± 6 3 

2695 160 ±32 6 

4850 50 ± 10 7 

4850 102 ±15 8 

4850 89 ± 10 9 

4850 58 ± 10 108 ± 5 36 ± 4 3 

10550 15 ± 4 60 ± 10 17 ± 5 3 



References: (1) WE NSS (iRengelink et al.lll997t). r2'l Texas (iDoudas et al.lll996D , (3) pre sent work, (4) NVSS 

(|Condon et al.llT998|). (5 VLandecker fc Caswelll (|198a) (6) iRcich et al.' (199(f). (7) MIT-GB (Langston et alJ|1990D, (8) 

87GB (ICxregorv fc Condonlll991h . (9) GB6 (fGreeorv et al...l99f).l . (10) VLSS (,Gohcn et al..,2007.;i 



Table 4 
Positions and spectra of compact sources 

Source 1 Source 2 Source 3 



RA(2000) 
Dec(2000) 

A 

B 


19'' 51"' 52.1" 19'' 53" 3.1" 19'' 52" 40.8" 

29° 26' 7" 29° 10' 36" 29° 35' 47" 

5.05 ± 0.09 4.69 ± 0.07 4.58 ± 0.11 

-0.92 ± 0.03 -0.72 ± 0.02 -0.83 ± 0.04 




log(S[mJy]) = A + B log(r/[MHz]) 



Table 5 
Integrated polarization properties 



Frequency [MHz] 


Polarized Intensity [Jy] 


Fractional Polarization [%] 


1420 
2695 
4850 
10550 


0.44 ± 0.02 
0.59 ± 0.05 
0.29 ± 0.03 
0.18 ± 0.02 




12 ± 2 
26 ± 4 
24 ± 4 
24 ± 6 



Fractional polarization is computed using integrated flux densities corrected for compact sources and the superimposed 

HII region. See Tabled 
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Table 6 
Breaks in the spectra of pulsar wind nebulae 



SNR 


Break Frequency 


injected/cooling" 


reference 


Crab Nebula 


40 keV & 1000 A 


i 


Woltier et al. (1997) 


Crab Nebula 


14000 GHz 


c 


Strom & Greidanus (1992) 


W44 


8000 GHz 


c 


Petre et al. (2002) 


VelaX 


100 GHz 


c 


Wciler & Panagia (1980) 


G29.7-0.3 


55 GHz 


i 


Bock & Gaensler (2005) 


3C58 


50 GHz 


i 


Salter et al. (1989): Green (1994) 


G21.5-0.9 


- 30 GHz 


? 


Salter et al. (1989) 


G16.7+0.1 


26 GHz 


i 


Bock & Gaensler (2005) 


CTB87 


10 GHz 


c 


Morsi & Reich (1987) 


G27.8+0.6 


5-10 GHz 


c 


Reich et al. (1984) 


G106.3+2.7 


4.5 GHz 


c 


Kothes et al. (2006a) 


DA 495 


1.3 GHz 


c 


this paper 



i denotes spectral break due to a break in the injected electron spectrum; 

cooling 



denotes spectral break due to synchrotron 



Table 7 

Present day and birth properties of the six pulsars for which a reliable age can be found. 

individual parameters see section 4.2. 



For a description of 



SNR 




Pulsar 


T [yr] 


E [erg/s] 


To [yr] 


Eo [erg/s] 


t[yr] 


Crab Nebula 




B0531+21 


1270 


4.4 X 10^8 


320 


1.0 X 10"*" 


950 


3C58 




J0205+6449 


5370 


2.7 X 10^7 


4550 


3.8 X 1037 


820 


Gil. 2-0.3 




J1811-1925 


23300 


6.4 X 103f5 


21680 


7.4 X 1036 


1620 


Vela X (/3 = 2 


.0) 


B0833-45 


11300 


6.9 X 10^6 


110 


7.3 X lO^*' 


11200 


Vela X (/3 = 6 


.0) 


B0833-45 


11300 


6.9 X 10^6 


3000 


2.0 X 10-*" 


8300 


G106.3+2.7 




J2229+6114 


10460 


2.2 X 10^7 


6560 


5.6 X 1037 


3900 


W44 




B1853+01 


20300 


4.3 X 10^5 


13300 


1.0 X 1036 


7000 



Table 8 
Properties of DA 495 extrapolated from birth properties of the six pulsars in Table [T] 



Basis 




iDA495 [yr] 


TDA495 [yr] 


Srcq [mG] 


^tot [erg] 


Smax [mG] 


Crab Nebula (/3 = 


:2.0) 


100900 


101200 


0.45 


1.0 X 10^" 


0.98 


Crab Nebula (/3 = 


:2.3) 


47600 


47900 


0.74 


7.8 X 10^9 


0.86 


3C58 




84100 


88700 


0.51 


5.2 X lO^* 


0.22 


Gil. 2-0.3 




165000 


186000 


0.32 


4.5 X 10^8 


0.22 


Vela X (/3 = 2.0) 




93900 


94000 


0.47 


2.5 X 10^" 


1.56 


Vela X (/3 == 6.0) 




20000 


23000 


1.32 


3.8 X 10^" 


1.91 


G106.3+2.7 




149000 


155000 


0.35 


1.1 X 10^3 


0.32 


W44 




28800 


42100 


1.03 


2.3 X 10^7 


0.05 
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29° 45 




19" 53" 



19" 52" 



19" 51 



19" 53" 



RIGHT ASCENSION (J2000) RIGHT ASGENSION (J2000) 

Fig. 1. — Total intensity maps of DA 495 at five radio frequencies, and the corresponding IRAS GOfim image. The position of the neutron 
star is shown in each panel by the white star. Contour levels are at: 10 K (12cr) to 100 K in steps of 10 K at 408 MHz (1 K = 8.5 mjy/beam); 
0.5 K (9(t) to 5.0 K in steps of 0.5 K at 1420 MHz (1 K = 8.5 mJy/beam); 0.1 K (5cr) to 0.9 K in steps of 0.1 K at 2695 MHz (1 K = 
400 mJy/beam); 5 mJy/beam {20cr) to 55 mJy/beam in steps of 5 mJy/beam at 4850 MHz (1 mJy/beam = 2.4 mK); 5 mJy/beam (4(t) to 
40 mJy/beam in steps of 5 mJy/beam at 10550 MHz (1 mJy/beam = 0.5 mK); 20 MJy/sr to 28 MJy/sr in steps of 1 MJy/sr at 60fini. 



DA 495 - an aging PWN 



17 



29° 45 




19" 52" 

RIGHT ASCENSION (J2000) 



19" 51' 



Fig. 2. — Separation of the total power emission at 1420 MHz into a non-thermal component (a = —0.85 — shown as black contours and 
shading) and a thermal component (a = —0.1 - white contours). Black contours arc from 30 mjy/beam (lOcr) to 170 mjy/beam in steps of 
20 mjy/beam. White contours arc from, 10 mjy/bcam (Scr) to 40 mjy/beam in steps of 5 mjy/bcam. The angular resolution is 2f45 (see 
Section 3.2). 
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Fig. 3. — Integrated radio continuum spectrum of DA 495. Flux densities taken from the literature are displayed in gray and those determined 
in this publication in black. The left panel shows a spectrum constructed from the measured flux densities of Table |2] The indicated line, 
corresponding to a = —0.59, is a fit to all plotted points. The right panel shows the spectrum after allowance for the contribution of compact 
sources and the superimposed H II region. The solid lines corresponding to a = —0.45 and a = —0.87 are fits to points below and above 
1 GHz respectively (see Section 3.2). The dashed line corresponds to q = —0.59, the fit from the left panel. The lengths of the arrows 
denoting the upper limits are not intended to indicate error estimates. 
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19" 53 



9^ 51' 



RIGHT ASCENSION (J2000) RIGHT ASGENSION (J2000) 

Fig. 4. — Polarized intensity of the emission from DA 495 and its surroundings at four frequencies, shown by contours and grayscale. The 
orientation of the E-field is shown by the superimposed vectors. The location of the pulsar is indicated in each panel by the star. Contour 
levels are at: 1.5 mjy/beam (4(t) to 12 mjy/beam at 1420 MHz; 6 mjy/beam (2(t) to 42 mjy/beam in steps of 6 mjy/beam at 2695 MHz; 
1.2 mjy/beam (60-) to 9.6 mjy/beam in steps of 1.2 mjy/beam at 4850 MHz; 0.8 mjy/beam (2a-) to 5.6 mjy/beam in steps of 0.8 mjy/beam 
at 10550 MHz. 
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RIGHT ASCENSION (J2000) 



RIGHT ASCENSION (J2000) 



Fig. 5. — Upper panels: Maps of rotation measure between 4850 MHz and 10550 MHz, at the resolution of the 2695 MHz map, 4.'3 (left) 
and at a resolution of 2.' 45 (right). Blaek eontours represent rotation measure, with values 140, 180, 220, 260, 300, 340, 380, and 420 rad/m^. 
White eontours indicate polarized intensity at 10550 MHz at the appropriate resolution. Lower panels: Left: Map of rotation measure between 
the four frequency bands of the CGPS at 2'45 (left). Contour levels are the same as in the other RM maps. Right: Map of polarized intensity 
at 4850 MHz shown by contours and grayscale, with levels 1.2 mjy/beam to 9.6 mjy/beam in steps of 1.2 mjy/beam. Overlaid white vectors 
indicate the intrinsic magnetic field directions calculated using the RM map at the top right. 
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Fig. 6. — Map of total intensity at 1420 MHz made by merging the CGPS 1420 MHz map with a map taken from the NVSS (see Section 
4.2). The resolution is about 45". Contour levels are at 0.8 (Scr), 1.6, 2.4, 3.2, 4.0, 4.8, and 5.6 K. Overlaid white contours indicate the 
polarized intensity at 4850 MHz with levels at 7, 10, 13, 16, and 19 mjy/beam. 
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Fig. 7. — Foreground rotation measure for all pulsars within 15° of DA 495 as a function of distance from us in kpc. See Section 4.4.1 for 
source of data. 
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Fig. 8. — Sketch of a possible magnetic field configuration inside DA 495, which contains a dipolc and a toroidal component. Dashed 
contours represent the polarized emission as observed at 4850 MHz. 
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Fig. 9. — Sketch of the dipolc field configuration adopted as a model of one component of the magnetic field in DA 495. The field lines have 
been assumed to have a dipole configuration that is stretched vertically by a factor of 2.5. The direction of the magnetic field IB is indicated. 
The spin axis of the pulsar should be aligned with the dipole axis. The viewing angle 9 is defined as shown. In Section 4.4.4 we consider this 
configuration viewed from the right. 
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Fig. 10. — The required electron density Ue as a function of viewing angle 8. The dashed lines in both panels correspond to the line 
of sight through the center of the dipole (Fig. ^, the lower line for the observed rotation measure of +170 rad/m^ and the upper line for 
+200 rad/m^, which is the observed value corrected for the foreground RM. The solid lines represent the line of sight through the dipole at a 
distance of 1.2 pc from the center, looking through the top of the dipole (left) and the bottom of it (right) as displayed in Fig. |9] The three 
lines were calculated from the top for rotation measures of +350, +300, and +250 rad/m^. The small dotted rectangles delineate the range 
of parameter values for which the model can provide a reasonable match to observed values of RM. 



O 
Q. 







1 





o 

Q. 






2 1 0-1-2 

Width [pc] 



2 1 0-1-2 

Depth [pc] 



Fig. 11. — Orientation of the dipole field inside DA 495, as determined in Sections 4.3 and 4.4. On the left the projection of the dipole field 
onto the plane of the sky is shown. On the right the line of sight orientation is indicated with view from the right hand side. 



